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ABSTRACT 


Since Newton, the stumbling block in celestial mechan- 
ics has been the three-body problem. Only restricted cases 
have yielded solutions. This paper describes a device, the 
"Lie-Series, "that first appeared in Lie’s work on analytical 
transformations; Grdbner has shown that they can be used 
to solve systems of differential equations by applying dif- 
ferential operators to known functions or to invert systems 
of analytical functions. The series are applied to Kepler’s 
problem of an undisturbed planet round the sun (two-body 
problem), to the study of perturbations, and to the process 
of obtaining the characteristics for any general dynamical 
problem. 
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ON THE APPLICATION OF LIE-SERIES 
TO THE PROBLEMS OF CELESTIAL MECHANICS* 


by 

Kaxl Stumpff 

Goddard Space Flight Center 


DEFINITION AND PROPERTIES OF LIE-SERIES 

The main problem of celestial mechanics is the study of the motions in space of n mass 
points that attract one another according to Newton’s law of gravitation. This problem can be 
solved for n - 2, but for n > 3 no general solutions exist that are valid for any values of the 
masses, initial conditions, and range of time. Since Newton, many attempts have been made to 
construct solutions for special cases, for instance, the solar system or the sun, earth, and moon 
system. The difficulty of solving the n -body problem justifies introducing any new mathematical 
method that gives deeper insight and facilitates solution. 

Recently, a new method has been proposed by Wolfgang Grobner (Reference 1). He uses the 
"Lie-Series" in solving a restricted class of differential equations. It is not yet clear that this 
new method is better than others; it has not been sufficiently used. But it is worthwhile to give a 
summary of the technique. The Lie-Series were originated by the celebrated Norwegian mathema- 
tician Sophus Lie, but they appear only incidentally in his works on analytical transformations. 
Their usefulness in solving problems like those of celestial mechanics was ascertained by Profes- 
sor Grobner. 

Suppose that 


*,<*) = e. (z,, z, z n ) 

are n functions of the n complex variables z k , and are holomorphic (i.e., can be developed into 
regular and convergent power series) within a certain neighborhood G around the position 

C = (5i- Q ■ 


* I his paper originates from lectures given at the Goddard Space Flight Center, Laboratory for Theoretical Studies, in February and 
March, 1967. 
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Define the operator D as 


d = *vo bt ^ 2 ( z > af; + ••• + e „(o-§r ■ (!) 

12 n 

This operator can be applied to any function f(z) that is holomorphic in G. Applying this operator 
repeatedly on f( z ) gives 


D 2 f = D(Df ), D 3 f = D(D 2 f) , 


and so on. Define a Lie-Series as 

1 t t D + -jrD 2 t •••] f( z ) . (2) 

v-0 


and write it symbolically 


L(z, t ) 


US* 


f(z) 


L = e tD f(z) 


This series is holomorphic in the region G and converges for all values of 1 1 1 smaller than a 
certain positive number T. The proof is left to the mathematicians. We shall assume that there 
is always a positive radius of convergence, no matter how small. 

The fundamental properties of the Lie-Series can be derived from the algebraic behavior of 
the operator D . It is enough to write down the principal aspects of the algebra. The following rules 
can easily be demonstrated: 


D [f j ( z ) t f , ( z )] 
D [c f ( z )] 

Dc 

D[f,(z) ■ f 2 (z)] 
D"[f, (z) • f , ( z )] 


Df ! (z) T Df 2 (z) , 
c Df ( z ) , 

0 . 

f, Df 2 (z) + f 2 Df j (z) , 


II 



a = 0 


f 2 (z) . 


From these, rules for the Lie-Series can be generated, for instance: 


et [c,f,(z) + 


2 n f n(z)] 


, t D 


f l( z ) 


D 


f ( z ) 

n v 7 


e tD [f ,( z) • f 2 (z) • • • f n (z>] 


[e tD f ,(z)] • [e tD f 2 (z)] • • ‘ [e tD f„(z)] 
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These two relations illustrate the following remarkable theorem: 


"The Lie-Series of a sum or a product of different functions is equal to the sum or the 
product of the Lie-Series of the components." 

This statement generalizes to the important rule 

e tD P(f v f a f„) = P(e tD f 1 .e tD f 2 , , e tD f n ) , 

i.e. the Lie-Series of a polynomial P is equal to the polynomial of the Lie-Series. This law is also 
valid for any analytical function of the z k . If we let z stand for z lt z 2 , ...» z n and F(z) for 
f( Z j, z 2 , ..., Z n ), we may write 

e tD F(z) - F(e tD z) or e tD F[f(z)] = F[e tD f(z)] ; (3) 

i.e., the symbol e tD and the function symbol F can be interchanged. This "rule of interchanging" 
is one of the most useful features of Lie-Series. The proof follows from the assumed uniform 
convergence of the Lie-Series within the region of their validity. Consider the simple example 
where z is a single complex variable (the operator D = d/dz). Then 

Dz - 1, D n z = 0 for n - 2, 3, . . . 


and 


e tD z - (1 + tD)z - z + t 


The rule of interchanging gives 


F(e tD z) - e tD F(z) 


or 


F(z + t) - £ td F<V) ( z ) 

v=0 


i.e., the well known Taylor series, which may be considered as a special Lie-Series. It is easy to 
extend these considerations to functions of several variables. 


The above simple application of the rule of interchanging shows that Lie-Series can demonstrate 
certain analytical properties in a very direct manner. This is still true wljen solving systems of 
first-order differential equations. 


Consider the equations 

dz . 

i 

dt 


0 iO) ; 


i = 1, 2 n 

(z) = (z 1P z 



2- • ’ Z n) 
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The functions 0.( z ) are holomorphic in a certain region G, and do not explicitly contain the 
independent variable t. Such a system is called "autonomous." (Non-autonomous systems can 
easily be reduced to autonomous ones, so that restriction to autonomous systems is not essential.) 
To solve the system in Equation (4) let the differential operator D be 

D = + SaCO^*'-- + ^ n (Olf7 


where 


(O = («!• 5 2 . ... £ n ) 

is a position in the region G that provisionally may be considered as variable but independent 
of the time. Then 


Z; - 4>i (£> t) - e tD 5; (i - 1, 2, , n) 

are n holomorphic functions of time t and of position i which, for t - 0, take the values 


Differentiating z. with respect to time gives 


= e-D Si 


according to the rule of interchanging. But because 


D£i = «.(0 


therefore, 


ITT = e tD ^(0 = e.(e‘ D s) = ^(z) 


i.e. the functions 


( 5 ) 


are solutions of Equations 4, with the initial conditions z.(0) = £. . The £ i are to be considered 
as variables during the application of the operator D, but after having performed all partial 



differentiations connected with this application they must be replaced by the constant values 
prescribed by the initial conditions of the problem. 

If the system of differential equations is not autonomous, then 

dzj 

~dt = 0 i( z .t) - 

where the functions 0. are holomorphic in a certain region g of the complex variables z t , z 2 , . . . , 
z n , t. Now introduce, only by a change of notation, the new variable z 0 = t , which is the solution 
of the equation dz 0 /dt = 1 , with the initial condition (z 0 ) 0 =o . This allows the nonautonomous 

system of n equations to be replaced by the autonomous system of n+1 equations 

dz i 

dt = 9 i{ z v) . (i = 0, 1, 2, .... n) 

with 0 O ( z) = l. This system can be solved using the differential operator 

3 V -1 3 

D = + 4, 0i(O 



in the same manner as before. 

The Lie-Series is applied here to a simple problem. Let A be the time, or another independent 
variable replacing the time, and let z(A) be a real or complex function of A. Then, if a 2 is a 
positive constant, 


z" + a. 2 z - 0, 



( 6 ) 


is the differential equation of a harmonic oscillation with the period 2 n/a. This problem can be 
solved directly by classical methods; however, we shall here use Lie-Series. The second-order 
differential equation is replaced by a system of two first-order differen tial equations: 

dz _ __ 

dA = y ~ e \ ( z - y) - 
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We let (£, 77) be a certain position in the (z, y)-space, and use the operator 



D 


77 aj -^ 2 C 


3t7 


Solving Equation 6, with the initial conditions 


z(0) 



for x = 0 , gives 




= ^1 + XD + AjfD 2 + •••J 


It is easy to apply the operator on £, as in this case D has a simple linear form. Indeed, we get 

DC = V 

D 2 C ~ D77 - ~CL 2 C , 

D 3 C- = -a 2 DC = ~a 2 V , ... 


and generally 


D 2n £ = (~l) n a 2n C 


D 2n + 1 C = (~l) n a 2n 7 / 


Therefore, the solution z is given by 


Z 


x 2 2 

C + X77 - 91 a £ 


3! 


2 

a^rj 


x 4 , 

+ 41 a 4 C 


(■' 


\ 2 a 2 

2! 


4! 


~ • • •) + 77 (x - 


3! 


X 5 a< 

5! 




or 


Z 


£ cos aX + sin 


(7) 


with the initial conditions' 


z(0) = C , 



( 8 ) 


This is the well-known general solution of Equation 6, 

ijX-n. 

6 

Y ■ ■ 



Suppose that the general solution of a system of differential equations 


dz . 

*dt~ = e i ( z > 


Si A 


D 




and 




the ^ being the coordinates of a certain point of the region G, in which the functions 0.( z) are 
holomorphic. Then every analytical function F(z) of the solutions z k is holomorphic in G, and 
may be written as a Lie-Series 


F(z) = e tD F(S) 


22 [f(o] 

= n 


This means that we can formally write down the development of any function F(z) without having 
explicitly calculated the solutions z k . Among these functions F, there are some for which 

DF(C) = o . 

These are called "characteristics” or "integrals" of the problem. Indeed, if F is a characteristic, 
then 


F(z) = e tD F(0 = F(C) 


i.e., F( z ) is constant on every trajectory z = z(t) which solves the problem. The last section of 
this paper discusses methods of constructing all the characteristics of a given system of differential 
equations. 

APPLICATION TO KEPLER’S PROBLEM 

The preceding section presented symbolically the general solution of a system of first-order 
differential equations with the aid of Lie-Series. Suppose that the equations are of the form 

dz i 

dt = 0 iO) - 0 i( Z l • Z 2 ’ ’ Z n) > 

7 


1, 2, ... , n 






the z i being complex or real coordinates, and the d i being functions that are holomorphic within 
a neighborhood G around a certain position (£) = M. The solutions of this system 

can be expressed by the Lie-Series 

(t'to)D 

z i = e C t , 


where 

D = 0 l(O 3^7 + + 0 n(Z) 


is an operator, and z. (t 0 ) = ^ are the initial conditions. All operations refer to £ as a variable, 
but afterwards i must be replaced by its constant value £=z(t 0 ). 

We may use this method to solve Kepler's problem of the undisturbed motion of a planet 
around the sun. The heliocentric rectangular coordinates of the planet being x 15 x 2 , x 3 , the 
differential equations of the problem are 

X. +— j - 0 , r = j/S X . 2 , i = 1, 2, 3, (9) 


if the units are canonic, i.e., if they are determined so that the gravitational constant k = l, and 
the sum of the masses Mtm = 1. These equations may be written as a system of the first-order 
equations: 

X . ” 

l 

u . = 

1 



If the initial conditions for t = t 0 are 


N 




( u 0 


Vi 


the differential operator will be 

D = 




3 £i_ _J_ 


P -- 


and the solutions which correspond to the conditions are 
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Carrying out the partial differentiations involved in the first of these formulas (the second is not 
absolutely needed) gives the well-known development of the undisturbed planetary coordinates in 
a power series of intermediate time r = t - t 0 . Unfortunately, this development is useful only 
for small intermediate times. For large r , higher-order terms in the series are needed, the 
expressions for the coefficients of r n become complicated for large n, and the convergence of 
the series is poor. However, if we introduce independent variable X. , where 


dA _ JL 
dr r ’ 


(ID 


the operator D becomes linear. The resulting solutions converge for all values of \ and may be 
written in closed form. Indeed, with pseudo-time A as independent parameter (x' = dx/dA) , 
Equations 9 change into linear ones, x ; " + a 2 x. = /?. = const. , a 2 = l/a = const. , a being the 
major semiaxis of a conic section; and, instead of Equation 10, 


- -„2 


a 2 x . + /3. 
1 ^ 1 


Another interesting method (Reference 2) of solving this problem makes use of the Levi-Civita 
transformation of the plane (x,,x^ -coordinates, 

X 1 - ^ 1 ~ ^>2 ’ X 2 ~ ^ 1 ^ 2 ’ ' 


which leads to the equations 


e, + 


= 0 


4a 


const . 


that are identical with Equation 6, and have the general solutions given in Equation 7, 


Vi 

E. - C- cos oA + — sinaA 
=* 1 ® l co 


with the initial conditions of Equation 8: 

£ 4 <0) = Ej , ^|(0) = Vi 
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Kepler's problem of the undisturbed motion of two bodies in a heliocentric system requires 
six integrals. The integral constants are the six elements of the orbit 

a, e, T; i, fl, co . 

The first three of these are independent of the choice of coordinate system; the last three 
are not. The first three are essential, as they determine the shape and the spatial extension of 
the orbit and its relation to the current of time. The last three, however, are merely directions 
for imbedding the orbit into an arbitrary given system of reference. Therefore to solve Kepler's 
problem it is sufficient to reduce the system of differential equations to one of the third order, 
using only variables that are independent of the system of coordinates. For instance, we may 
solve equation given on page 201 of Reference 3: 


for the radius vector r, provided that for r = t-t 0 = o the initial conditions r(0), f(0), ?( 0) are 
given. When the pseudo-time A is introduced, the differential equation takes the linear form 


l-r" 1 

a - ” - - const. 


which has the well-known solution 


r - a + b cos a A. + c sin aA 


( 12 ) 


a, b, c being constants depending on r 0 , r 0 , r 0 . To obtain this solution by using Lie-Series, 
replace the one third-order differential equation by first-order equations. For instance, by the 
system 


r ' 

- <p{ r, 

4 >. '/<) 

<t>‘ 

= Wr. 

<f>, 4>) 

r 

= -a 2 

<t> ■ 


i — 

r r , 


which is accomplished using Equation 11 (which connects A and time). Then the differential operator 


3r 0 + * ~° 2 4 > 30 0 


3 


D 


+ 



where r 0 , <p 0 , <p 0 , r 0 are the values of r, <p, «/>, r for \=o, gives the general solutions 


r 


e^ D r, 


4> - e^ B <P 0 , 


+ = , 


r = e XD r 0 


It is sufficient here to consider the first and the last of these formulas. Considering the first, we 
have 

r = ^1 + \D + ~2T D 2 + • • r Q , 

hr 0 - 4> o — r o > 

D 2 r 0 = TH> 0 = 0 O = r" , 

D 3r o = rV'o = - a2< t>0 = _a2r o • 

and generally, for n > l , 

D 2n r Q - (-a)"' 1 !-", D 2 "' 1 r„ = (- a )" -1 r' . ( 13 ) 


Therefore, 


A 2 . 

A 3 „ . 

A . 4 2 „ 




1 + TT r o -~3T a T 0 

- 4T a 2 r 0 + 




( A 3 

. A 5 . 


'X 2 

A 4 , 

A 6 4 \ 

^ "IT 1 +_ 5T a 

“ + r o( 

,2! 

“ 4 1 o- 2 

+ ~6T 0,4 - •••) 


or 


r 


r o 



sin 



(1 


cos a X) 


which is Equation 12 in another form. If we introduce the c -functions 


sin a A. 

c 0 - cos a\ , c, - - 


1 — cos a. A. _ cl \ ~ s in a. A. 

(a \) 2 ’ ° 3 = (aX ) 3 


we may write the solution for r in the elegant form 

r = r 0 + r 0 C 1 K + r 0 C 2 K2 



Jf/. 
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In a si m ilar way, we derive the connection between time r and pseudo-time X. by solving the 
differential equation 


dr 

d\ “ r 


by the Lie-Series 


T 


e t 


o 


As r Q = 0, 


= r o - d2t o = Dr o D»t 0 = Dn ' lr o ■ 

Using Equation 13 for the operation of D n on r gives 

r = r O X + r 0 C 2 X2 + r 0 C 3^ 3 ' 

the "main equation" of Kepler's problem, which is equivalent to Kepler's equation 

n(t-T) - E-esinE 

Trying to solve Kepler's equation has fascinated generations of astronomers and mathemati- 
cians. Countless methods have been used, especially to overcome the difficulties when orbit 
eccentricity e approaches 1 (the parabolic case.) Today, interest in this problem has notably 
decreased because the main equation (the regularized form of Kepler's equation) avoids these 
difficulties and is applicable to orbits of all eccentricities without singularities. 

Solving either Kepler's equation or the main equation may be considered an example of 
inverting analytical functions. We have equations of the form 

0 = t - T - f(E) (Kepler’s equation) , 
t = t - t Q = g(^) (Main equation) , 

and the problem is to find the reverse relations 

E = F(0) , \ = G(t) . 


12 



I 


1 


This inversion is possible under the following general conditions, 
functions 




Given a system of analytical 



Vi - ^i< x ) - ( x i » x 2- - X n) - 

where the functions <p { are holomorphic in the environment of a position 

Vi = <^(£) = ^(^l* ^ 2* * * * ’^n) ■ 

and the Jacobian matrix 



is different from zero at the position £, then the system of functions is reversible; i.e.j there 
exists the solution 


X. - <t.(y) - <j>.( yi> y 2 1 .... y n ) , 


with the inverse Jacobian matrix 


r 1 


(*i0 



The Lie-Series are helpful in solving this problem. Indeed, the x. as functions of the y L must 
satisfy the system of partial differential equations 

Bx . 

= ®iu( x ) 

with the initial conditions 


for y. 


V, 


and it can easily be shown that it is possible to write the solution of this problem in the form of 
Lie-Series. 

For simplicity and with special attention to the problem of solving Kepler's equation, we shall 
consider here only the case of one function 

y = f(x) ; y ~ V f° r x = £ . 


m 
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We can solve for the inverted function 


X - <£( y) 


provided that 3y/3x f 0 for x = . In this case, we can write 


dy 

= *( x ) 


f'(x) 


as a common differential equation, the non-singular inverse is 


dx _ 1 

dy - f ' ( x ) > 


provided f'(0 f o . 

Then by Lie-Series, with 


D 


l d 
f'CO d<f - 


x = e<y T >> D £ 


(14) 


is the solution of Equation 14 with the prescribed initial conditions. 
Apply this rule to Kepler's equation 

y = x-esinx 

Then y* = 1 - e cosx is different from zero for x = 0 and e < 1, and 


1 - e cos df d^ 


Applying this operator repeatedly to we obtain 

l 

^ ’ l-o cos £, ’ 

n 2 A = _ e sin ^ , 

u = (1 - e cos § ) 3 


D 3 <f 


e cos E, 


3 e 2 sin 2 E, 


(1-e cos f) 4 (1-e cos f) s 



Therefore, 


. i n 



x - € + (y ~v)°£ + ~21 (y~v') 2 D 2 ^ + yr(y — t?) 3 d 3 ^ + 


and, for <f = 77 = 0, 


1 - e 3! 


(1 ~e) 4 


This is indeed one of the simplest special series solving Kepler’s equation. 
It is interesting to apply this method to the main equation 


r - r 0 \ + r ' c 2 \ 2 + r" c 3 X 3 . 


If 


D = 


1 a 1 a 


r o 2K o 


r o ^0 


denotes the solving operator, then 


X = e TD X„ 


r r 2 

- 1 + tD + I D 2 + ••• X Q 


and 


^ = 0 


DX, 


D 2 X_ = D 


1 /_ r o' \ _ _ jV 

r ° \ r oV = ‘ r 0 3 


D 3 X„ 


-#) 


1 T o 
r 0 r » 


3r 0 2 


Therefore, 


x = 


r o 2! r 3 


r 3 / r o" 3r o 

3! I r ~ 5 

\ r o r o 
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or, introducing the dimensionless quantity 



and multiplying \ by r 0 /r gives 



When these considerations are extended to systems of several functions of several variables, 
the process of inversion is complicated as we are dealing with partial differential equations. 
Without going into details, note that this problem requires, not one operator D but n operators, 
such that 



i = 1 


and these operators must fulfill the conditions of reversibility 


which correspond to the conditions of compatibility, 


3 2 x. 
1 

3 yj 


3 2 x. 

i 

3y u 3yj 


between second-order partial derivatives. 


APPLICATION TO THE PROBLEM OF PERTURBATIONS 

Lie-Series are helpful in symbolically representing the solution of any system of first-order 
differ entiateequations. A theorem of Lagrange, applied to a more complicated problem, demon- 
strates the usefulness of Lie-Series. 

Suppose that three bodies (for instance, the sun and two planets) are moving according to 
Newton's law of gravitation. Let the masses be nq, m 3 , m s , and let q 15 q 3 , q s be the vectors 
leading from m 3 to m 5 , m s to m 1} m 1 to m 3 respectively, and 

Qj “ q 2 ’ ^3 _ ’ ^5 ^6 


be the corresponding velocity vectors. This system of relative motions requires 12 integrals of 
motion. Between the q's there exist the vector relations 



^3 + q s 


= o 


q 2 + + le 


Therefore, only four independent vectors need be solved for, e.g., , q 2 , q 3 , q 4 . Lagrange states 

that this three-body problem is solved if we know, as functions of time, nine independent functions 
of the q ; that are invariant under coordinate transformations. This solution being achieved, the 
remaining three integrals can be obtained directly or by simple quadratures. The problem of 
finding the nine first integrals is called the "reduced problem of three bodies" (see Reference 4, 
or Chapter XII of Reference 5). If 


pik = (p ; • ?k) = x ; x k + y^k + z ; z k ■ 

(the scalar products between any two of the six vectors q, inclusive of the cases i T k), it is 
sufficient to solve for the nine independent quantities 


Pll’ P 33’ PSS’ P 1 2 ’ P34’ PS6' P22’ P44’ P66 


(15) 


as functions of the time. These quantities may be called "fundamental invariants" of the problem, 
as they are independent of the system of coordinates and of one another, and form a symmetrically 
builtup system. These nine quantities are solutions of the differential equations 


_ 1 _ . 

2 Pi 1 


J_ . 

2 P 22 



I 

S’ 

1 

<-l 3 

n "3/2 

Pi 1 

Pll + m l p l ’ 

with ■< 

Pi 

= MjPn + 

^ 3 Pl3 

+ ^sPis 

+ m lP2 , 



- ^!P 21 + 

^3P 2 3 

+ M 5 P 2S 


(16) 


and six others are obtained by cyclic permutation of the indices 1, 3, 5,and 2, 4, 1 . 6 , respectively. 

The invariants Pj and p 2 (likewise p 3 and p 4 , p 5 , and p 6 ) on the right-hand sides can be expressed 
as functions of the fundamental invariants by a number of easily constructed relations between the 
p ;k (e.g., Pll + p 13 + p 15 " 0 , p 32 + p 34 + p 36 “ 0 , etc.). We could supplement this system of 
nine differential equations by those for all the other possible Pjk invariants. The total number 
of different invariants is 21 (for N bodies there are N • (2N+1) p’s). The system 


dp 


i k 


dt 


^ik(Pll’ P 1 2 ’ ’ P 56 ’ P66) 


e ik(P) 



may be solved by Lie-Series, if we use the operator 


D 


i e. 


3 


3p 


ik 


the sum containing 21 expressions. The solutions can be written 

p = e( t-t °) D 7T.. 

r ik ik ’ 

where ^ ik are the values of the p ik for t = t 0 . In discussing Lagrange's theorem it is enough to 
write down these developments for nine of the fundamental invariants (Equation 15) or for nine 
independent functions of these invariants. The procedure may be improved by introducing a 
pseudo -time s instead of t. 

This method of developing the solutions of the reduced three-body problem does not seem to 
have been attempted. Perhaps it is as intricate and cumbersome as any of the classical methods, 
but the short and elegant symbolic form of the solution may be more tractable for computer pro- 
gramming, as the right-hand sides of the differential equations are finite algebraic expressions of 
the invariants themselves. 

Another interesting, possible use for Lie-Series is in the case of slightly perturbed planetary 
orbits. Then, as in classical celestial mechanics, we may divide the motion of a planet into two 
parts: an intermediate orbit, (for instance, a Keplerian ellipse easily calculable with simple sub- 
routines) and a perturbation term proportional to the small perturbing mass. The differential 
operator that solves the problem may be written as the sum of two parts, 

D = Dj + D 2 . 

These two operators must fulfill some necessary conditions: 

1. Dj must be constituted so that the solutions 


are simple^andvwell known functions. This is always the case if D x represents the linear part of 
D, for instance if D x is the operator of a Kepler motion osculating the true orbit, and if a linearizing 
pseudo-time is used as independent variable. 

2. d 2 must be small compared with D x ; i.e., the coefficients 6 y of d 2 must be noticeably 
smaller than those of D x , at least within a certain range of time. If these conditions are satisfied, 



the problem can be solved by the Lie-Series 



When D=Dj + D 2 is applied to f , the f is considered as variable, and set equal to f (t 0 ) only after 
all differentiations have been performed. However, this exponential expression developed into a 
power series of r = t - 1 0 gives 

f = [i + ^( D i + » 2 ) + 4 ( D i + d 2 ) 2 + ■ ' f N ; 


and 


+ D 2 ^ 2 expands to Dj 2 + DjD 2 + D 2 D l + ’ etc - > 

because, generally, D,d 2 / d 2 Dj , i.e., the operators Dj and D 2 do not commute, except under very 
special conditions. Let us digress and consider briefly the commutation relations for the d's. 

It can be shown that in the reduced three-body problem D t and d 2 do not commute. Write the 
solutions of the three differential equations for Pll , p 12 , p 22 , using the differential operator 
D = D t + d 2 with 


2p>12 3 Pn 




a 

9 Pl 2 


2p l2 P l 1 


■ 3/2 


3p 


22 


D 


2 


Pi 


a a ~ 

a Pl2 +2p 2 ap 22 _ 


supposing that Pl , p 2 are functions of Pll , Pl2 , p 22 and known functions of the time (i.e., supposing 
that the motion of the perturbing mass m, is known). Then, for instance, 


lPll 2 P 1 2 > 

D 1 ( D 2 P l l) 

= D,(0) 

2 Pll = o - 

D 2 ( D lPll) = 

2 D 2 Pl 2 


therefore D t D 2 ^ D 2 D 2 . 

It is easy to derive the commutation conditions for any two operators d 2 and d 2 . Suppose there 
are two operators 
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Then 


D l®2 ^11 3zj (^21 3Zj + ® 


2 2 3z. 


+ ^12 3z 2 ^21 BZj + ^22 3z 2 + ) 


Be 


21 B 


+ e„ 


B 2 


Bzj Bzj 1 ^21 9z 2 


The product of the two operators consists of two terms, one containing first-order partials and the 
other containing only second-order partials. A simple calculation shows that, for the operator 


(D x . d 2 ) = DjD 2 - d 2Di , 


the sum of the second-order terms is zero. Therefore (d^ d 2 ) is a linear operator like D : and d 2 . 
This operator has the symmetric form 


(°i. = 22 ’ 

i = i 


where 


0 . 

1 




- e 


2k 



Dj and d 2 commute only if (D lt D 2 ), the "Jacobi bracket", is zero. The Jacobi brackets behave 
like the bracket expressions of Lagrange and Poisson, which are well known in celestial mechanics. 
The Jacobi brackets obey the rule 

((D,, D a ), D a ) + ((D 2 . D 3 ), Dj + ((D 3 , D*), D a ) = 0 , 


which is. also valid for the brackets of Lagrange and Poisson. 

Let uscreturn to the problem of perturbed motions. Let z = (z,,..., z n ) be an n -dimensional 
vector describing^the position of a system of bodies as a function of the time, and let i = , . . • , £ n ) 

be the position at t = 0 . Then 


- ‘( D . + °2) v 

- e 4, 


z(t) 


f(t, O 



will be the solution of the system of differential equations belonging to this problem, 
the exponential, we may write 


z(t) - 




r—l 

2^^ D r 2 d 2 dz 



Developmg 


v= 2 


CO 

J2 FT Dj V-a D 2 D a ” 1 


z + 


z = £ 


(17) 


On the right side z must be replaced by £ when all operations are performed. The first sum, 


f 0 c t . z ) = 2 Z^ tDiVz ~ e ' Diz 

v = 0 


contains all terms that depend on only and symbolize the intermediary orbit. The second term 
contains d 2 only once and at the end. The third sum has d 2 at the last position but one, and in 
certain terms at the last one too, and so on. This series is not in a convenient form, as the 
perturbations of different order are not separated. On the contrary, there are parts of the per- 
turbations of the first order in all sums except the first, and there are parts of the second-order 
perturbations in all sums but the first and the second. The perturbations of different orders can 
be obtained separately only by our totally changing the sequence of the terms; this is possible if 
the series converges absolutely. 

For simplicity, consider only the first-order perturbations. Select from each sum the terms 
that contain D 2 only once; i.e., from the sum beginning with v = a take the partial sum 


2ZS D i v ‘ aD 2 D r 1 2 

v- a 


Then the perturbation terms of the first order are completely collected in the double sum 


f i (t> z) = D r aD 2 D r 


l v + a + 1 


(v + a + 1 ) ! 

a ~ 0 v - 0 


D 1 a D 2 D 1 ‘ / z 


(18) 


Each term of this expression contains D 2 only once, but in every possible position amongith'el;. 
factors Dj 


2 — * — ’ — j *■ *■ — 

Similar considerations lead to the expressions f n (t, z) for the p er turbationsfafimgher 


order. 





Grobner proposes another method which gives the solutions 


'(t) - [ f 0 ( t - z ) + f 1(0 z ) + • • '] 


> = { 


in the form of an integral equation of the Volterra type. We can rewrite the complete partial sum 
of Equation 17: 


co 

L 


D 1 y -“D 2 D 11 ' 1 


CD 

L 


D, v D„ D a_1 z 


( v + a) ! 


by making use of the identity 


l v +a 


(v + a) ! 


■ f 


(t -T )^ 1 J 

(a - 1)! v\ dT 


which can easily be proved by repeated integration by parts, in the form 


r t co 

f V 1 (t T v 

J 0 Z_ (a-l)! vT D " D 2 D“-i z dr 


This is the partial sum with index a. Summing up all these expressions and changing the index a 
into a + l gives the total expression for the solution z(t) , 


z(t) 


f(t, z) 




-t °0 / CO 

I feSw 

a-0 \ v ='0 


dr 




Now we may write in the inner sum, 


co 



v = 0 


and, according to the rule of interchanging, 


’ 1 (o 2 D a ) 


<2 = 0 


d 2 d c 




D 2 D a z 


( z = f o) 


Therefore, 


L 

n =rv 


(t -T) 
n * 


a 


D 2 D a z 


V - ’ (t - r y 1 
D 2 TP D<l2 


D 2 f(t~r, z) 


and finally 


z(t ) 


foO- z > z = 4 



~ T , 


2) -f, 


dr 


This integral equation can be solved by iteration; this process will rapidly converge if D 2 is small 
enough. The first steps of the iteration will be 

tD , 

z o ' ^o( t ’ z ) z =C ' e z (z=C) 

r t 

Z 1 = f 0 + _T - Z )z = z n dT - 

Jo u 


and one can easily prove that the integral on the right side of the second equation is identical with 
the double sum in Equation 18 derived for the first-order perturbations. 


APPLICATION TO THE CONSTRUCTION OF THE CHARACTERISTICS OF 
DYNAMICAL PROCESSES 

Three previous sections applied Lie-Series to some problems of celestial mechanics: the 
two-body problem, the reduced three-body problem, and the special problem of computing general 
perturbations. 

Is this new method of solving dynamical problems an improvement? That depends on our point 
of view. The development of solutions in terms of power series is scarcely changed by the innova- 
tion. But the Lie-Series give short, pregnant solutions for any dynamical process— an obvious 
advantage. The Lie operators have a universal form that encompasses nearly all the types of 
linear operators used in former theories. The very interesting properties of the Lie-Series, 
especially the important rule of interchanging (Equation 3), allow the expressions for the solutions 
of dynamical problems to be reduced to simple and most lucid forms. This is helpful, especially 
in developing series for computer programming. 

This section discusses a method to obtain the integrals of motion of a dynamical process. 
Consider a concrete case: the reduced problem of two bodies. This involves three integrals— in 
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other words, three independent constant functions, which may or may not contain the time explicitly. 
In the reduced problems, the integrals are invariants with respect to coordinate transformations. 

If Xj , y , , and z, are the heliocentric coordinates of a planet, and x, = x,, y 2 = y 1} z 2 = z, 
the coordinates of the velocity, then the reduced problem of two bodies will be solved when three 
fundamental invariants are found (as functions of time): 

2 . 2 ■ 2 - 2 
Pn " x i T Vi T z i ' r ■ 

p 1 2 = x 3 x 2 + y j y 2 + z ) z 2 = rV cos • 

2 , 2 . 2 - \i2 

p 22 - x 2 - + y 2 ‘ + Z 2 - V 

where r is the distance from sun to planet, v is the magnitude of the relative velocity, and <p is 
the angle between the vectors of position and velocity. These invariants are the solutions of 
Equations 16 (with - p )I “ 3/2 , and m, - o): 

Pu = 2 p J2 . 

Pl2 

P 2 2 


P 2 2 P 1 l 


- 1 2 


o -3 2 

~P l 2 P 1 1 


The integrals of this problem are 

l . a 

— - 2p tl ' - p 2 , (the; integral of energy) , 
a( 1 _ e 2 ) r Pii P 2 2 — P i (the integral of areas) 

T = t-a 32 (E-e sin E) = f(t, a, e)(Kepler’s equation). 


These integrals provide three arbitrary constants a, e, and T. The first two integrals describe 
the geometrical form of the heliocentric orbit of the planet; the last integral describes the planet's 
motion in this orbit. The first two integrals show that certain functions of the variables p ik are 
true time-independent constants, denoted as "characteristics" of the problem. In the reduced two- 
body problem the characteristics are finite algebraic functions of the variables. The third integral 
behaves quite differently. This constant expression contains the time explicitly and involves the 
time-dependent function e, usually denoted as the "eccentric anomaly". This is a very simple 
case, showing that among the three integrals of the problem there are two characteristic functions. 


24 



Generally we can prove that among the n integrals of an n -order problem there are exactly n-i 
independent characteristics, i.e., n-l independent functions of the variables that are constant along 
each special solution of the problem and do not contain time explicitly. 

The proof that an n -order problem has n-l characteristics at most, is easy. Let 

D “ 0 i( z )^7T + ( z )^r~ 

1 n 

be the Lie-operator of the problem, and let F(z) be a characteristic function; then DF( Z ) ~ o. 

If F,( z ) , . . . ,F r (z) are r different characteristics, each analytic function g(f,, f 2 ,..., F r ) is a 
characteristic also, as 


DG 


5G 

1f7 df i 



0 


As each characteristic F.( z ) is a function of the n variables Zj , z 2 , ..., z of the problem, it is 
clear that no more than n independent characteristics exist. Indeed, if there were n^l functions 
f k of the n variables z _ , there would exist a relation g(f, , F 2 , ... , F n+1 ) = 0. Actually, there 
are no more than n-l independent characteristics. Assume that there are n characteristic 
functions 


w i - Fj(z) 


F (z) 

n v ' 


If these functions are independent, the Jacobian matrix 



must be regular. Therefore, it must be possible to invert the equations and obtain 

Z 1 ' ■ • ■ ■ ’ Z n " L,0) • 

But the functions (w ,, . . . ,w n ) are characteristics, as shown above, and therefore 


Dz. - d(z) - 0 

1 1 v ' 

i.e., the problem restricts itself to the trivial case 

z. S' ^.(z) = 0 , i - 1, 2, ...n 


which can be excluded. 
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Finally, we must prove that n-l independent characteristics really exist. Let us assume that 
the solutions of a system of n differential equations 

z . - r . ( z ) 

1 l v ' 


are 


z i = ( e ' Dz ,)z = - :: • 

Then, eliminating t from any two of these equations gives a characteristic. The inverse function 
of any analytic function 


y - /-■( x ) 


in the environment of a position r, = -4(z) , where <£'(x) / o. may be expressed by the Lie-Series 


X 


(y~r)D \ 
e X _ , 

i / *-s 


if 


D 


1 _d_ 

."(x) dx 


is an operator that is regular in the neighborhood of x = - . If we apply this rule to any two solu- 
tions of our problem 


z , = /.('• O ■ 
z u r M 4 - ■ ) • 


we obtain the reversal expressions 


t 





with the operators 


D 

l 


J_ d_ 

l. dt • 

* i 


-k 


_d 

dt 


( 19 ) 


which are regular in any region around t = t where the derivatives of ? i and i k do not vanish. The 
difference between these two expressions for t gives an equation 


.1 
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with constant values £ and independent of time. This formula represents a characteristic of the 
problem, as it is valid along each solution. Combining thus the first solution with the other n-l 
solutions gives n-l independent characteristics 



M z i- z K ; 0 ~ 0 


2 , 3 ,. 


( 20 ) 


Besides the n-i characteristics in Equation 20, there exists one independent relation (Equation 19) 
that gives t as a function of the particular z i% • e.g., 


t = 





This equation, t = f ( Zj , £ 15 . . . , £ n ) , plays the same role in the general problem as Kepler's 
equation does in the two-body problem. In this case, z 1 takes a favored position but can obviously 
be replaced by any other of the z n . 

The above device for constructing the characteristics of a dynamical problem does not quite 
satisfy our expectations, as each of the cases of Equation 20 contains only two of the n variables 
and more or less of the constants 5 ; . (Pertinent to this question there is an unpublished remark 
made in 1950 by K. Stumpff.) 

Let Cj • • • £ n = £ be the position of a vector z t . . . z n = z(t) at t = t 0 . Then, r = t - t 0 , 


“ z i ‘ z i T + 2 


Z . Z . 

2 _ 3 


+ 


where z . and its derivatives are taken at time t 0 . If z. is different from zero. 


^i( T - Z 1 Z n) 


as the derivatives Z; (n) are functions of the z themselves, and may be derived from the differential 
equations z. = 6>.(z). Reversing these equations gives, as before, n equations 

T = f. (z, {,) . 

Suppose that one of the (for instance £ n ) is zero. Then 

T = f n (z) 


is a function of the variables z only, 
gives n-l equations 


Substituting f 


for r the equations for f l , 





fi(z, Q 


f n( Z > 
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from which 


£i ~ £j(z) - const. , i = 1, 2, , n-1 , 

by inversion, provided that the position i is so chosen that in a certain surrounding region G the 
z x are not zero. These equations, which can easily be expressed by Lie-Series, are the charac- 
teristics in a convenient form, and they describe the geometrical feature of the trajectories. These 
n-l constants are supplemented by the main equation r = f n (z), which expresses the connection 
between the variables z and the time. 
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